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We give an improved estimate of primordial magnetic fields generated during cosmological first 
order phase transitions. We examine the charge distribution at the nucleated bubble wall and its 
dynamics. We consider instabilities on the bubble walls developing during the phase transition, 
ft is found that damping of these instabilities due to viscosity and heat conductivity caused by 
particle diffusion can be important in the QCD phase transition, but is probably negligible in the 
electroweak transition. We show how such instabilities together with the surface charge densities 
on bubble walls excite magnetic fields within a certain range of wavelengths. We discuss how these 
magnetic seed fields may be amplified by MHD effects in the turbulent fiuid. The strength and 
spectrum of the primordial magnetic field at the present time for the cases where this mechanism 
was operative during the electroweak or the QCD phase transition are estimated. On a 10 Mpc 
comoving scale, field strengths of the order 10~^^ G for electroweak and 10~^° G for QCD, could be 
attained for reasonable phase transition parameters. 
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I. INTRODUCTION 
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The origin of galactic and intra-cluster magnetic fields 
and the existence of a cosmological magnetic field are still 
unknown. It has been suggested that dynamo effects in 
turbulent fluids might amplify small seed fields exponen- 
tially 0]. However, it has been shown recently |^,^ that 
dynamo theory has difficulties explaining the galactic and 
cluster magnetic fields. Alternatively, these fields may 
be the result of compression of a primordial field which 
would then permeate the universe as a whole. In this sce- 
nario, the required strength of the primordial field on a 
Mpc scale would be much higher, -Bprim ^ 10~^^ — 10~^ G 
at present Q , compared to the case where the fields have 
been amplified by a dynamo, i?sccd ^ 10^^"' G. 

It is thus of some interest to investigate mechanisms 
which could produce primordial fields and estimate their 
strength and spectrum. Scenarios involving inflation can 
lead to large coherence lengths, but the predicted ampli- 
tude is in general very small (5|j^. Magnetic seed fields 
can be created by the slightly different dynamical re- 
sponse of the negatively and positively charged compo- 
nents of a quasineutral fluid such as in the Biermann bat- 



tery mechanism 0. This thermoelectric effect operates 
as shocks form with the nucleation of bubbles during the 
QCD transition but the created fields are also quite 
small. More recently, magnetic seed fields caused by the 
turbulent flow near the walls of bubbles which are be- 
ing nucleated and grow during a cosmological first order 
phase transition have been estimated for the QCD |^ and 
the electroweak (EW) transition. These two publi- 
cations make different assumptions about the nature of 
the charge layers and the fluid flow at the bubble walls. 
As a result, the seed field strength estimates in Ref. 
are larger by some 14 orders of magnitude compared to 
the ones in ref. |0 when each generated field is com- 
pared to the appropriate equipartition strength for each 
transition. 

Scenarios where magnetic fields are created directly 
from the Higgs fields involved in the EW transition have 
been discussed in Refs. [l^Jl^. Estimates of magnetic 
fields generated during phase equilibration of a complex 
Higgs field caused by bubble collisions in a first order 
phase transition have been given in Ref. JTst . 

In this paper, we estimate the strength of primordial 
magnetic fields generated during a cosmological first or- 
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der phase transition. We rederive the charge density dis- 
tribution at the bubble wall in Sec. II. We then discuss 
instabilities in the bubble growth as well as their damp- 
ing due to a finite shear viscosity and heat conductivity 
caused by particle diffusion in Sec. III. In Sec. IV, we 
show how the fluid flows associated with the bubble wall 
instabilities together with the charge density distribution 
at the wall will lead to the generation of magnetic seed 
fields. In this section we determine the wavelength range 
in which the magnetic fields are initially excited and the 
amplitude of these seed fields. In Sec. V, we discuss the 
transition to the non-linear regime and how the fluid tur- 
bulence can amplify the initial seed fields. We apply our 
scenario to the EW and the QCD transition in Sec. VI 
and estimate the resulting strength and spectrum of the 
extragalactic magnetic field at present time. A summary 
of our findings is contained in Sec. VII. We will use nat- 
ural units and the CGS system throughout the paper. 



II. SURFACE CHARGE DENSITIES AT THE 
QCD AND THE EW PHASE BOUNDARIES 

An unmagnetized plasma may develop magnetic fields 
if the spontaneous separation of electric charges give rise 
to a net current. Separation of charges will occur during 
a first-order EW transition through the development of 
net baryon number or net top-quark number gradients 
in or at the EW phase boundaries, regardless of whether 
baryogenesis does or does not occur at the EW transition. 
During a first order QCD transition net baryon number 
discontinuities of appreciable magnitude may develop at 
the phase boundaries leading to the spontaneous separa- 
tion of charges. This charge separation has been previ- 
ously considered by Ref. |H| for QCD phase boundaries 
and by Ref. for EW phase boundaries. 

The first order phase transition proceeds from a high- 
temperature ("h") to a low-temperature ("1") phase. For 
the EW transition this would correspond to the phase 
with a vanishing and non-vanishing vacuum expectation 
value of the Higgs field, and for the QCD transition to the 
quark-gluon phase and the hadronic phase, respectively. 

In the QCD transition the baryon number susceptibil- 
ity in the high-temperature phase is larger than that in 
the low-temperature phase such that in chemical equi- 
librium there is a net baryon number density contrast, 
Rb = {rib) eg /{nl) eg, whcrc {n^)eg is the net baryon 
number for the quark-gluon phase in equilibrium (high- 
temperature) and {n.l)eg is the same for the hadron phase 
(low-temperature). Ri, ^ 10—10'^ depending on the tran- 
sition temperature. However, chemical equilibrium can 
only obtain for small bubble wall velocities, vj, <^ 1, and, 
in general, the net baryon number contrast can exceed 
Rb significantly. 

The width of the baryon number excess layer and the 
baryon number density contrast between the quark-gluon 
phase and the hadron phase can be estimated by studying 



the steady-state solutions for the net baryon density in 
the wall rest frame. Requiring that the net baryon flux 
into a thin shell in the quark-gluon phase at the phase 
boundary be zero, we find that 

V"6 RbJ Lb 

The first term in Eq. (|l|) represents the net flux of baryon 
number across the phase boundary where v'j^ is the ther- 
mal nucleon speed in the hadron phase and T,h is the 
probability for a nucleon approaching the phase bound- 
ary from the hadron phase to dissociate into three quarks 
and pass over into the quark-gluon phase. The second 
and third terms in Eq.(Q) give the diffusive baryon flux 
in the quark-gluon phase away from the phase bound- 
ary and the advective ffux towards the phase boundary, 
respectively. Here Iq is the quark mean free path and 
Lb is a characteristic length scale of variation in baryon 
density. For large bubble velocities, Vb ~ 0.1, and/or 
small baryon number penetrability, S/j <c; 1 jl^ , the 
flux across the phase boundary is negligible and chem- 
ical equilibrium between the phases cannot obtain. In 
this case, an exponentially decaying baryon number layer 
with width Lb — Iq/vb develops at the phase boundary 
in the quark-gluon phase while far away from the wall 
baryon number is homogeneously distributed. The out- 
of-equilibrium density discontinuity at the phase bound- 
ary can then be much larger than its equilibrium value, 
("b/^-fe) ^ Lib- In the opposite limit, e.g. Wb <C 1 and/or 
Yih — 1, equilibrium is approximately maintained such 
that net baryon density is roughly homogeneous in both 
phases. Nevertheless, there still exists a baryon number 
discontinuity at the phase boundary with density con- 
trast {nb/n\) ~ Rb- 

For the case of a first order EW transition, let us flrst 
consider the case where baryogenesis occurs at energy 
scales higher than the EW breaking scale. In this case 
a cosmic baryon number asymmetry (implying a top- 
quark number asymmetry) already exists in the high- 
temperature phase through, for example, the existence of 
an asymmetry in (baryon - lepton) number which is con- 
served by baryon number violating EW sphaleron pro- 
cesses. Then very similar arguments to net baryon num- 
ber in the QCD transition hold for net top-quark number 
in the EW transition. Top quark number densities are 
thermodynamically suppressed in the low-temperature 
phase due to their finite rest mass such that in equilib- 
rium {n^/nl) — Rt ~ 2. By writing a similar equation 
to Eq.(|l|) for top quarks at the EW phase boundaries, 
and with vf^^ ~ 1, 'Stop — 1, and Rt only slightly larger 
than unity, one finds that only for relativistic wall veloc- 
ities, Vb — 1, chemical equilibrium cannot be maintained. 
Since we consider iif, <C 1 (see below), net top quark num- 
ber will be in equilibrium which implies that top quark 
number density will vary by a factor of order unity over 
the extension of the EW bubble wall. 

If EW baryogenesis is operative during a first order EW 
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transition, the cosmic baryon asymmetry will be gener- 
ated either in front of the wall in non-adiabatic (thin 
wall) baryogenesis scenarios or within the wall in adia- 
batic (thick wall) baryogenesis scenarios jl^. Note that 
for the simplest extensions of the standard model of EW 
interactions one expects thick walls, L^aU ^ 1/?" [0- 
In this case, a baryon asymmetry (carried by u, d, s, 
c, and b quarks) will smoothly vary from zero in the 
high-temperature phase to its cosmic average value deep 
within the wall where baryon number violating processes 
shut off. 

In the QCD transition, the net baryon number is asso- 
ciated with a net positive charge density. In particular, 
Pc = (2/3)en„— {l/3)end in the quark-gluon phase in the 
limit where the strange quark mass is large compared to 
the temperature, and pc ~ eup in the hadron phase where 
Uu, nd, and Hp arc net number densities for up-quarks, 
down-quarks, and protons, respectively. Globally this 
net positive charge is neutralized by a negative charge 
density due to the net electron number. However, locally 
on the scale of the baryon discontinuity/inhomogeneity 
at the phase boundary electrons are free streaming, e.g. 
le ~ [(aVT2)(g/7r2)T3]-i > U, such that the thin 
excess baryon number layer/discontinuity at the phase 
boundary cannot remain neutral. In the absence of elec- 
tric fields, riu — Ud (for n„), while weak 
equilibrium in the presence of a positive electric poten- 
tial requires Ud > «-« and leads to screening of the charge 
density by the conversion of positively charged up-quarks 
into negatively charged down-quarks. 

Consider the grand potential density for up- and down- 
quarks in the presence of an electric field, i?. 
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where A = (47re^T^)~-'^/^ in this model. We note here 
that the actual Debye screening length in a multi-species 
relativistic plasma is \d — {i-Ke^gcT^)"^^^ where gc is 
the statistical weight of relativistic charged particles with 
mean free path I ^ max[A£), Lh]. 

The Debye equation can be solved for the expected 
steady state net baryon number distribution in the QCD 
transition and the net top-quark/baryon number distri- 
bution in the EW transition. For the QCD phase bound- 
aries we find a surface charge density 
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(5) 



(2) 



Here the first term arises from a charge density on the 
wall (the width of the wall is of order ~ 1 fm <C A/? , ap- 
proximately a QCD-color screening length) and is spread 
over approximately one electric Debye screening length, 
whereas the second term results from a much smaller 
charge density on the scale of the baryon number inho- 
mogeneity in the quark-gluon plasma, if, 3> A^i. The 
positive charge at the phase boundary in the quark- 
gluon plasma is balanced by a negative charge density 
at the wall in the hadron phase, resulting in a thin elec- 
tric dipole layer. The baryon density at the wall in the 
quark-gluon phase, ^ n^, is approximately ~ 2nf, 
halfway into the transition when chemical equilibrium 
is maintained, but can become appreciable towards the 
end of the transition, ^ Rbfib- Here fif, is the cos- 
mic average baryon density, fii, ~ 40ry5((7/100)r^ with 
775 ~ 5 X 10^^^ the cosmic average baryon-to-entropy 
ratio. 

In the EW transition there are two separate cases de- 
pending on the width of the EW phase boundary, in 
particular either thick walls Lwaii ^ Ad, or thin walls 
iwaii ^ A_D. The resulting surface charge density, if EW 
baryogenesis is not operative, is 



where (f) is the electrostatic potential, pi and qi are chem- 
ical potential and electric charge of species i, and Q is 
given to lowest order in p[ = pi — qicf). Baryon number 
density rif, and charge density pc can be obtained from the 
net quark number densities Ui = —(dQ/dpi)T.(i>, which 
yield 



nb = ^^^(a^L 



Pc^-eT^{2p',~p'^) 



(3) 



Further demanding weak equilibrium in the case where 
Pe = Pv — ^ (electrons and neutrinos are free streaming) 
gives Pd = Pu + g4>. The Poisson equation together with 
Eq. (|^) can then be used to derive the Debye screening 
equation 



(4) 



en. 



'■top 



1 



en. 



top 

h 

''top 



L 



wall ■ 



1 I Ad , if iwaii ^ Ad , 



top 



if -^waii ^ Ad 
(6) 



and is seen to be suppressed by (AD/iwaii)^ for contin- 



uous thick walls. Here n'^^^ ~ (l/6)nf, > 



ni^p for all but 



the largest wall velocities, Vb — ^- The magnitude and 
spatial extension of charge densities when EW baryoge- 
nesis does occur is of the same order as that given in 
Eq.(i. 

To summarize, we parametrize the charge density at 
the bubble wall by 



(7) 



which extends over a length scale 1^ ^ fc/Tc around 
the waU. In the QCD case 77 - 5 x 10"^" - 10-^ and 
/c ~ 1 for the dominant electric charge density resulting 
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from the baryon number discontinuity right at the phase 
boundary, and for the EW transition, 77 ^ 10~^^ and 
fc ~ 40, where we assume thick walls, Lwaii ~ 40/T. 



III. INSTABILITIES OF BUBBLE GROWTH AND 
VISCOUS DAMPING IN THE LINEAR REGIME 

During cosmological first order phase transitions, in- 
stabilities may develop as bubbles grow. Hydrodynamic 
instabilities can occur when the transport of latent heat 
is dominated by the fluid flow. These have been inves- 
tigated for the QCD transition and the EW tran- 
sition in the small velocities limit, for cosmological 
detonation fronts p9| ] , and for general first order transi- 
tions in the limit of very small or very large velocities [po[ . 
These bubble wall instabilities may be damped by finite 
viscosity and heat conductivity due to the diffusion of ra- 
diation on small length scales. To account for this damp- 
ing we will use the approach in Refs. |21 2^ for length 
scales larger than the radiation mean free path. 

We note that, depending on the parameters of the 
phase transition, the heat transport may be dominated 
by diffusion instead of convection. In that case, bubbles 
larger than the radiation mean free path would likely be- 
come unstable to "dendritic growth" We discuss 
below under which conditions this may happen. 

The high and low temperature phases for either the 
EW or the QCD transition can be described by the fol- 
lowing two equations of state for pressure Pi (T) , enthalpy 
density Wi{T), energy density pilT) {i — h,l), and tem- 
perature T: 



pi{T) = [wi{T) + L]/4, 
Ph{T)=Wh{T)/A, 



pi{T)^[iwi{T)~L]/A, 
Ph{T)^?,Wh{T)/A, (8) 



where m,(T) = {2t:'^ / 'i'o)giT'^ {i = h,l), and 5, is the 
number of relativistic degrees of freedom in the plasma. 
Furthermore, S = {wh — wi)/wh — L/wh > 0, with 
Wi = Wi{Tc) {i = h,l) and L the latent heat. The crit- 
ical temperature Tc is defined by pi{Tc) — ph(Tc) or, 
equivalently, by equating the free energy densities in the 
two phases. We assume that charged and strongly inter- 
acting particles have approximately the same velocity v 
and can thus be described in the one-fluid approximation. 
The fluid carries a conserved quantum number, namely 
baryon number, with a flux given by — nU^. Here, n 
is the proper baryon number density and C/'' ~ (Ijv) is 
the four velocity (we assume non-relativistic flows). 

In order to obtain the fluid flow and electromagnetic 
(EM) flelds during bubble nucleation in a first order 
phase transition, one would have to solve the following 
equations in the presence of a phase boundary simulta- 
neously: Energy-momentum conservation, T'"' -i, = 0, 
where Tt"" = {p + p)^^ - pgt^" + r^'' + T^^ is the 
energy-momentum tensor including non-ideal contribu- 
tions r'^" from a finite viscosity and heat conductivity 
caused by the diffusion of photons and neutrinos |pi| as 



well as EM contributions Tg^; baryon number conserva- 
tion, = 0; and, finally. Maxwell's equations. 

A full solution of this set of non-linear equations 
would require extensive numerical calculations. Analyt- 
ical approximations are usually obtained by linearizing 
the equations around their zeroth order equilibrium so- 
lution. This has been done in the literature for various 
limiting cases: For example, in Ref. |2^, the case of only 
one phase was considered and dispersion relations and 
damping rates for MHD waves were obtained. In this 
work. Maxwell's equations were treated in the ideal MHD 
approximation (assuming infinite conductivity which is 
an excellent approximation for the evolution of the mag- 
netic field B in the plasma under consideration ^), i.e. 
the hydromagnetic equation 



atB = V X (v X B) 



(9) 



was used. 

For an ideal fluid (i.e. vanishing viscosity and heat 
conductivity, r^" = 0) and no EM fields, the remaining 
linearized equations for the perturbations in pressure, p' , 
and velocity, v', are [|17[|1S] 



{dt + v-V)p' + -V -V = (10) 

{dt+v- V) v' + -Vp' - Iv V • v' = , 

w 3 

and have to be solved in both phases. Both the zeroth 
order solutions for a planar bubble wall, represented by 
the unperturbed quantities v and p (or w) and the per- 
turbations v' and p' have to be matched at the phase 
boundary by requiring continuity of energy and momen- 
tum flow and the transverse velocity across the interface. 
The latter condition comes from assuming equilibrium in 
the presence of a finite shear viscosity. It has been shown 
in Ref. that growing instabilities seem to be possible 
only if the bubble wall velocity Vh satisfies 



Vb < Vc 



21 1/2 



(11) 



where Th is the temperature in the high-temperature 
phase. The treatment in Refs. 0,11^ is only adequate 
in the limit Vb <C Wcrit- It was also pointed out pC| ] 
that the QCD transition is a borderline case since Vb — 
Vciit — 0.03, whereas the EW transition is unlikely to 
fulfill Eq. ( pT| ) in the initial stage of fast bubble growth. 
However, since Vb cx {T^ — T^)/T^ @, the condition 
Eq. (|ll|) may eventually be satisfied if — > Tc- Indeed, 
Ref. argued that the EW bubble walls may slow down 
significantly (by 1-2 orders of magnitude) if reheating 
is significant. For reasonable phase transition parame- 
ters, Vb < 10~^ during this stage, and Eq. ( pT| ) can be 
met. It is therefore possible that hydrodynamic instabil- 
ities develop both in the QCD transition and during the 
late stages of the EW transition for which the approach 
in Refs. [ pTlp^ ] can be applied. In the following, we will 
assume that this is indeed the case. 
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Here, we are interested in solutions of the general lin- 
earized equations for two phases including viscosity and 
heat conductivity and EM fields. This is still a quite 
complicated problem and one usually considers the sim- 
pler limiting cases of either one phase of a viscous fluid 
or two phases of an ideal fluid. In the general case, we 
expect a certain wavelength range over which there are 
growing instability modes, i.e. the growing instability 
dominates over the diffusive damping. Since we are only 
interested in an order of magnitude estimate of the re- 
sulting initial seed fields and their coherence scale, we 
can restrict ourselves to the case where the damping rate 
is small compared to the instabilitygrowth rate. We can 
then describe these modes by Eq. (|o|) and verify a poste- 
riori that the damping rate 7d of the corresponding bub- 
ble wall instability is small compared to its growth rate 
7inst- Seed fields can then be generated by these growing 
bubble wall instabilities. We will verify that backreac- 
tion effects of these seed fields onto the fluid flow in the 
growing instabilities are negligible also. 

Let Vh and vi be the modulus of the unperturbed veloc- 
ity of the high- and low-temperature phase, respectively, 
in the bubble wall rest frame. We assume that the flrst 
order transition proceeds as a weak, non-relativistic de- 
flagration, i.e. Vh.vi < Cs, and (5 ~ (w; - Vh)/vh < 1, 
where Cs — l/\/3 is the speed of sound, and get [p7|p8| 



7mst - Svhk/2 . 



(12) 



This holds true in the range of wavenumbers fcmin ^ k ^ 

fcmax, with 



2/RS, 
Swhvl/a, 



(13) 



Here, the first condition comes from requiring that 7d be 
larger than the expansion rate v^/ R of a spherical bub- 
ble of radius i?, whereas the second marks the bound- 
ary to the regime where the bubble wall is stabilized 
by the surface tension, cr ~ T^. The hydrodynamic 
instability therefore can become operative for bubbles 
with a radius R large enough such that fcmin < fc,„ax. 
This is certainly the case for the fully developed bubbles 
near percolation which we consider in this paper. For 
these bubbles, i? ~ i?. 



perc 



fbTH, where rn ^ 2.3 x 
10*'(5/10)-i/2(r/100MeV)-2cm is the Hubble radius, 
and g is the number of relativistic degrees of freedom in 
the plasma. For the QCD case, /& ~ IQ-^ - 10"^ ||,|6|, 
and for the EW case ~ IQ-^ - 10"'' 

The growth rate in Eq. ( p^ has to be compared with 
the damping rate 



9t . 
7d ~ — mm 

55 



k Tj., 

Tr 



(14) 



where gt is the statistical weight for particles efficiently 
transporting momentum, either neutrinos {gt = 5.5) or 
photons [gt — 2), and = ((cr,.) gT^ /tt"^)^^ is the mean 



free path of these particles, and {ar) is the correspond- 
ing average cross section. For fcr,. < 1, Eq. (^4|) was 
obtained from Ref. |^2|, and for kry > 1, the damping 
rate is approximately independent of k and proportional 
to 1/Tr ||2^. Associated with the damping rate of insta- 
bilities is the Reynolds number of the fluid flow at length 
scale 1/fc, 



R{k) 



(15) 



where vf{k) is the typical fluid velocity at that length 
scale. In the linear regime of perturbations, Vf{k) ^ v^. 

For the EW transition, the dominant processes 
governing diffusion of radiation turn out to be EM 
bremsstrahlung, e + I —>■ e + I + j and photon-electron 
pair production, e + ^ —^ e + l~^l~ , where I stands for any 
charged lepton in the plasma. These processes are influ- 
enced by plasma effects. With a the fine structure con- 
stant, Uc the density of charged particles in the plasma, 
and mp ~ {Anaric/STy^^ ~ 0.5T{g/100y^^ the plasma 
mass of the electron at temperature T ^ 1 MeV [g8[ , the 
relevant cross sections can be written as (ar) = (ctem) — 
(SaVm^) ln(6.3r/mp) ~ 3 x 10-^{g/100)-^/T^. This 
leads to the Reynolds number for EM viscosity, 



i?EM(fc) =i4x 10^ 



xfbSvf{k) 



VlOO/ 
k ■ 



lOOGcV 



(16) 



where mpi is the Planck mass. Using the numerical pa- 
rameters fb^& few xlO~^, 6 ^ 10~^, and ~ — 
vi ^ 0.01 in the percolation stage of the EW transi- 
tion il,||,|9|, we get i?EM(fcmi„) > 10^. 

For the QCD transition, the damping rate Eq. 
is dominated by neutrino diffusion for k < (temTu)' 
where tem and Tj, are the mean free paths of EM radi- 
ation and neutrinos, respectively. Using (ar) = (cr^) ~ 
2.1G|T^, where Gp is Fermi's constant, we obtain for 
the Reynolds number for neutrino viscosity 



RAk) 



ioM4 



1/2 / T, 
10/ V 100 MeV 



fbSvf{k) ■ 



(17) 



Using the parameters /{, ^ 10 , S ^ 10 , and Vh ~ 
0.1 |29] for the QCD transition, we get i?;y (fcmin) ^ 1 and 
-REM(fcmin ) ~ 10^^ 

Combining Eqs. (Oh and (H), we obtain 



7d 

7inst 



0.2 

Svh 



and instabilities are expected to be damped in 
the wavenumber range 5((7/10)(gt/5.5)~ii?(fcinin) ^ 
k/krnin < 0.2{g/W)-\gt/5.5)R{k^in)/{Svh)^. From this 
we see immediately, that, at least in case of the EW 
transition, damping is unlikely to be important at the 
largest length scales over which instabilities can develop. 
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Damping probably only plays a role on length 
scales which are a few orders of magnitude smaller than 
1/fcmin, in a window of relatively small logarithmic width 
of ~ 2 log [{gt / 5. 5) /Svhig/ 10)] - 1.4. In contrast, for 
the QCD transition, damping could become important 

around kniin, if Ru{kmin) ^ 1- 

We note that since 7inst(fcmin) — Vh/R is roughly the 
convection rate, 7d(fcmin)/7inst(fcmin) ^ 1 corresponds to 
the case where transport of latent heat is dominated by 
diffusion instead of hydrodynamic flow. Therefore, de- 
pending on the exact phase transition parameters enter- 
ing Eqs. (p^ , ([l7|) , and (p^, this situation could occur in 
the QCD transition, as remarked in Ref . , but not for 
the EW transition. 



IV. SEED FIELDS FROM HYDRODYNAMIC 
INSTABILITIES 

Seed fields associated with instabilities will be gener- 
ated when positive and negative charges are displaced 
relative to each other and, therefore, the one fluid ap- 
proximation cannot be used. However, the resulting EM 
currents, j ~ pcV ^ erjT^v, correspond to very small 
conduction velocities, Vc, of negative relative to positive 
charge carriers, ~ <SC v. Therefore, the subsequent 
evolution of the seed fields due to dynamo effects can 
very well be described within the MHD approximation. 
Note that seed fields are necessary since the MHD equa- 
tion @ is homogeneous in B. We now estimate the seed 
fields generated by the hydrodynamic instabilities. 

Let the planar, undisturbed phase boundary be located 
at z = with the low-temperature phase at z < and 
the high-temperature phase at z > 0. We can choose 
the instability to have a wavevector in the x direction 
and the perturbation of the wall can be described by 
Zuj{x,t) = Zoexp(7inst^ + ikx). This corresponds to a 
nonpropagating instability with Im7i„st = 0. To lowest 
non-trivial order in Vh and 6, and for k <^ fcmax, the 
solutions of Eq. ( ]To| ) take a relatively simple form: 

Vft. = Svh (e^ - ie^) kz^{x, t)e^^^ 12 , 

Pft. ^ -Swhvlkzyj{x, t)e-''y2 , (19) 



Svh 



^e*W2 + + ^e,)efc- fcz^(a;,i)/2. 



p'l = -5whVj^kzw{x, t)e '^72 , 

where e^; and are unit vectors in the x and z direction, 
respectively. Note that Eq. ( [l9| ) only holds for distances 
from the wall which arc much larger than the perturba- 
tion amplitude of the wall location, i.e. z ^ Zw{x, t). It is 
seen from Eq. (19) that the instability perturbs the fluid 
flow in the low-temperature phase up to a characteristic 
distance 1/Sk, whereas in the high-temperature phase the 
fluid flow is only perturbed up to a distance 1/fc. For the 
wavelength regime under consideration, Eq. (p3|), these 
scales are typically much larger than the width ^ fc/Tc 
of the dipole layer. 



Eq. (|l^ illustrates that, for a given mode with k = 
fce^;, the X component of the perturbed velocity field has 
a discontinuity at z — 0: 



(vJi - V;) • e.j; = -i6vhkz^{x, t) . 



(20) 



The dipole layer with charge density given by Eq. (0) 
thus leads to a net current in the x direction. 



J, ~ eriT^8vhkz^{x,t)e-^^\'\'f^ , 



(21) 



provided fc/Tc ^ Zw{x^t). We can now superpose all in- 
stability modes from Eq. ( |l9[ ) with k in the z = plane 
and write the wall displacement spectrum as (z^) = 
^[zw{k)Y'dh\k. Using Ampere's law without the dis- 
placement current, V x B = 47rj, and Eq. (|2l|), after 
performing a three dimensional Fourier transformation, 
we arrive at 



B,{k) ^ AeriT^Svhfc (fcminfc)'/' Km 



(22) 



for the seed field spectrum (-Bg) = J[i3s(fc)]^o?ln fc. 
Strictly speaking, Eq. ( p2| ) is only reliable for Zu,(fc) < 

fclTc. 



The instabilities Eq. (19) grow non- linear when the 
perturbed velocity spectrum 



«/(fc)~(5l/V(fcmi„fc)'/'k»(fc)| 



(23) 



satisfies w/(fc) ^ w^, where, analogous to the seed 
field spectrum, the flow velocity spectrum is written as 

^Wj^ = /[t!/(fc)]^dlnfc. Furthermore, the seed held B^ 

and its source currents j give rise to an additional force 
acting onto the fluid of the form j x B. This term is of 
second order in the perturbation amplitude and was ne- 
glected in Eq. (|l|). We now compare its magnitude with 
the leading force term Vp' which is of first order in the 
perturbation. By using Eq s. (p9| ), Ampere's law and an 
equation analogous to Eq. (|23D for |Vp'|, we obtain 



\jBs\ {enfcf n 
IVp'l " 



Vh Wh 



S'/\f{k)., 



(24) 



on a length scale 1/fc. This is small compared to unity 
in the linear regime where Vf{k) ^ Vh- In addition, the 
electric field Es{k) induced by the growing magnetic seed 
field, Es{k) ^ [dBs{k)/dt]/k ~ 6vhBs(k), is smafl com- 
pared to the zeroth order field E ^ erjfcT^ caused by the 
electric charge layer, Es{k)/E ~ Vf{k) <C 1. We there- 
fore conclude that in most cases backreaction effects of 
the seed fields onto the medium are negligible in the lin- 
ear regime. Therefore, a conservative estimate of the seed 
field Bs is given by assuming z^(fc) ^ fc/Tc', 



Bs{k) ^ ierjfcT^S^^\f{k) > Ae7jT^Sv,j! 



T2 



1/2 



(25) 
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Being less conservative, one can extrapolate the first ex- 
pression in Eq. ( ^ ) up to the transition to the non-linear 
regime, where Vf{k) ~ Vh, and obtain the optimistic esti- 
mate Bsik) ~ erjT^S^^^Vhfc- In the next section we will 
show that once the non-linear regime is reached, enough 
turbulence in the fluid is expected to build up and the 
seed fields can be amplified by dynamo effects. 

For comparison, the thermoelectric effect discussed in 
Ref. ||] for the QCD transition leads to seed fields of the 
order ~ O.OIw/Tc/i?pcrc which is smaller than Eq. (25) 
for k 3> Aimin. This can be seen by rewriting Eq. (25) 
as B,{k) > ivvhfciTc/Rperc)ik/k^i^y^\ Other mech- 
anisms have been investigated in the literature starting 
from thermal field fluctuations on a length scale ~ I/Tc, 
Bff^{k ~ Tc) ~ 8ttT^. On larger scales, k < Tc, the mag- 
netic field spectrum was argued to behave as a power-law, 
Bth{k) ~ (87r)i/2T2(fc/rc)P. Several options have been 
discussed for p: Based on stochastic arguments applied 
to the order parameter involved in t he p hase transition, 
Ref. ||n| argued for p = 1, and Ref. (f|] for p = 1/2. In 
contrast, treating the resulting magnetic dipole moments 
as stochastic variables leads to p = 3/2 In any case, 
by comparing with Eq. ( ^5|) we see that for p ^ 1, very 
roughly, on a scale k ~ fcmin given by Eq. (^3|) , our mech- 
anism leads to comparable or stronger fields already in 
the linear regime. 



V. AMPLIFICATION OF SEED FIELDS BY MHD 
EFFECTS 

The seed fields discussed in the previous section can be 
amplified by exchange of energy with the turbulent fluid 
flow. This non-linear MHD regime is usually investigated 
by numerical simulations |3^ . The magnetic field energy 
Em typically grows up to equipartition with the fluid 
flow where it takes the value 



f>2 

rjii 2 



(26) 



with Beq being the corresponding equipartition field. The 
fluid velocity in the turbulent regime can be estimated by 
assuming that the latent heat release L ~ Swh is at least 
partially converted into turbulent motion. This results 
in Vf ~ (5^/2 which is probably somewhat larger than 
Vf '^Vh in the linear regime. 

We note in passing that the resulting strong equipar- 
tition fields Beq on bubble scales do not cause a problem 
with recent bounds from big-bang nucleosynthesis 
This is because fields on scales smaller than the comov- 
ing horizon size at neutrino-decoupling are expected to 
be significantly damped by neutrino diffusion by the time 
nucleosynthesis commences p2| . 

The required amount of amplification for the seed field 
in Eq. (E5|) is given by the inverse of the ratio 



Bsik) 

Ben 



-1/2 



>10 



-17 



too GeV ) /f,(5i/2 



Vh 



(27) 



which is > 10-24(7^/10-14) for the EW and > 
10-24(7^/10"^) for the QCD transition, for the param- 
eters used above. Note, that in case of the EW tran- 
sition, BJBeq is much larger than what was estimated 
in Ref. 110] if it is scaled to the same -q. The reason 
is that our instability analysis revealed the existence of 
a monopole layer of an EM current on the bubble wall, 
whereas in Ref. ||l^ a dipole layer was assumed. 

The central question in the theory of magnetic field 
amplification in a turbulent fluid concerns the time scale 
over which energy equipartition is achieved and the spec- 
trum of the resulting magnetic field. Here we use the 
analytic approach and physical arguments developed in 
Ref. 1^ which was recently applied to the problem of 
galactic field generation |Q. Once the bubbles start 
to collide, vorticity is expected to develop which, on 
a scale k is of the order kv'{k) and corresponds to a 
turnover rate of equal size. This turnover rate becomes 
larger than the instability growth rate Eq. (12) once 
v'{k) ^ Svh- If the Reynolds number R{k) ^ 1, one 
therefore expects the fluid flow to turn turbulent be- 
fore the instabilities become non-linear, which occurs for 
v'{k) ~ Vfi- At that point, as long as the velocities 
involved are smalled compared to the speed of sound 
Cg ~ l/\/3, one expects a Kolmogorov type of spec- 
trum to develop. The typical velocity on a scale k is 
then given by w/(fc) ~ Vf {ke/kY/^ Q. This holds for 
kc k fcvis — kcR{kc)^/^ ^ where kc is the wavenum- 
ber for which the eddy turnover rate ~ kcVf is equal to 
the inverse of the percolation time tporc — Rpcrc/vh, giv- 
ing kc ~ Vh/RpcicVf ~ l/i?pcrc. The viscous scale /cvis, 
at which dissipation becomes important, is basically the 
scale of the smallest eddy in the turbulent flow. 

Neglecting resistive damping for the moment, the 
growth rate of the magnetic field energy during the MHD 
regime is given by [S| 



dE 



M 



dt 



27a/ - / Vfik)dk, 



(28) 



(29) 



with 



where the velocity spectrum is again written as \V'j ^ 

J[vf{k)]'^d\nk. Thus, 7m is basically the maximum 
turnover rate of the eddies in the flow. For a Kolmogorov 
spectrum, 



2/3 

7Af ~ Vfkc ( ^ ) ~ Vfkc [R{kc 



11/2 



(30) 



The resistive scale kr is defined as the wavenumber 
where the resistive damping rate ~ ~fc^/(47rcrcond) 
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caused by the finite conductivity (Tcond — lOTc p^ ] 
becomes comparable to the growth rate 27m- Thus, 
kr ~ (47rcrcond7M)^/^ ^ Kis[RM/R{kc)]^/^, whcrc Rai = 
47r(27r/fcc)w/crcond ^ 10^ (5/100)"^ i?EM(fcc) is the mag- 
netic Reynolds number. The magnetic field spectrum 
cuts off exponentially for k > kr- It has been shown in 
Ref. [^4| that 7m is constant and Em grows exponen- 
tially until 2jmEm becomes comparable with the tur- 
bulent power P ~ gT^fvj{k)dk ~ gT^k^vj In R{kc) 
which is also the dissipation rate due to viscosity. At 
that point, the Alfven velocity va ^ {2EM/gT^Y/'^ be- 
comes comparable to the turbulent velocity at the vis- 
cous scale, ^/(fcvis) — VfR{kc)~^^^, which imphes Em — 
E'^R{kc)~^^'^ ■ Since the magnetic field on scales k with 
kvA ^ 7m — fcvisW/(fcvis) cannot grow, the magnetic field 
energy should be concentrated around k^is by the time 
when P ~ 2-fMEM- From then on, 2^mEm — P is 
roughly constant, so that Em only grows linearly in time 
until it reaches E'^j. Therefore, at the same time, jm 
has to decrease until it reaches ~ v/fcc- By comparing 
with the second equality in Eq. (|30|), this means that the 
high wavenumber cutoff in the velocity power spectrum 
decreases until it reaches fee- The viscous length scale 
grows due to the extra viscosity caused by the energy 
drain from the smallest eddy to the magnetic field. The 
same happens to the magnetic field spectrum, so that in 
the end all the magnetic energy E'^j is concentrated at 
k = kc, i.e. at the largest wavelengths. The time scale 
over which this process occurs can be estimated as 



dE 



M 



peg 



1 



2jmEm P kcVf 



pcrc ■ 



(31) 



We have thus verified that amplification up to equiparti- 
tion is possible within the percolation time. Furthermore, 
the magnetic field is expected to be coherent on the scale 
kc, and we can write B{l/kc,Tc) ~ {8TrgSy/'^T^ . 



VI. RELIC EXTRAGALACTIC MAGNETIC 
FIELDS 

On length scales r ^ 1/kc we will again assume a 
power-law behavior cx {kcr)~P for the magnetic field. 
Since kc is a macroscopic scale compared to the corre- 
lation length of the order parameter of the phase tran- 
sition which is ^ 1/Tc, the large scale fields should be 
determined by the randomly oriented magnetic dipole 
moments, leading to p ~ 3/2 [^l]j9|Jl0t] . For the magnetic 
field on a comoving scale r at redshift z this eventually 
leads to 



S(r,z)^10--(l-Hz)V,^/V/2(4) 



Tc 



100 GcV 



6 

-3/2 



lOMpc 



1/4 

-3/2 



G. (32) 



We could have followed previous authors by choosing 
a fiducial value for r of 1 Mpc, which corresponds to the 



scale that collapsed to form a galaxy. The estimated field 
would then be appropriate for seeding the galactic fields. 
However, it was recently suggested that magnetic fields 
on scales smaller than the Silk scale ~ 10 Mpc are likely 
to be damped by photon diffusion at recombination [ p2[ . 
Thus, Eq. (^H) probably only applies for r > 10 Mpc 
and leads to an extragalactic magnetic field spectrum 
B{r,z = 0) ~ 10-29 (r/10Mpc)-3/2G for the EW tran- 
sition, and B{r,z = 0) - 10-20(r/10 Mpc)-^/^ q foj. ^j^g 
QCD transition for the parameters used above. In the 
latter case our mechanism thus predicts a field strength 
on a 10 Mpc scale which is roughly comparable to the field 
expected from the Biermann battery mechanism acting 
on intergalactic distance scales during large-scale struc- 
ture formation [^ . 

For comparison, the stochastic models mentioned in 
Sec. Ill lead to a field strength B{r, z) ~ {%ttY/'^T^{1 + 
z)'^{rTQ)~P, where Tq ~ 2.7 K is the present temperature 
of the cosmic microwave background. For both the EW 
and the QCD transition this gives BjlO Mpc, z 
10-19 G and - 10-32 G for p 1/2 || and p = 
respectively. 

In Refs. |37 3I] it was pointed out that hydromagnetic 



I'M 



turbulence in the early universe might shift a given mag- 
netic field power spectrum near equipartition with the 
fluid to length scales larger by a factor / which can be 
as large as a few orders of magnitude. In that case, the 
length scale r in the above estimates should be substi- 
tuted by r//. At large scales characterized by a given r, 
this could increase the fields by a few orders of magni- 
tude. 



VII. CONCLUSIONS 



We have discussed the generation of magnetic fields 
during the growth of weak deflagration bubbles nucle- 
ated in a cosmological first order phase transition in some 
detail. This proceeds essentially in two steps: First, in- 
stabilities can arise which can be described by linear per- 
turbation theory. Together with a finite charge density 
along the bubble walls these instabilities cause EM cur- 
rents and thus magnetic seed fields. While currently we 
are unable to decide unambiguously whether such insta- 
bilities form, for reasonable phase transition parameters 
this is at least plausible. Taking into account damping 
due to a finite viscosity and heat conductivity caused by 
radiation diffusion we showed that the resulting bubble 
instabilities can grow on length scales somewhat smaller 
than the bubble radius at percolation. Once these per- 
turbations grow non-linear, the fiuid is expected to turn 
turbulent and the seed fields can be amplified by MHD 
effects in the form of a dynamo. Qualitative physical 
arguments show that equipartition of the magnetic field 
energy with the kinetic energy in the turbulent motion 
can be achieved and the field spectrum can be concen- 
trated at length scales not much smaller than the bubble 
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radius at percolation. Random superposition of the mag- 
netic dipole moments associated with the bubbles leads 
to large scale fields which can be estimated and have a 
characteristic spectrum. On a lOMpc comoving scale, 
field strengths of the order 10~^^ G from the EW tran- 
sition and 10^^" G from the QCD transition could be 
attained at z = for reasonable phase transition param- 
eters. Finally, hydromagnetic turbulence after the phase 
transition may further enhance the large scale fields. 
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